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Abstract 

The proton energy spectrum and the angular distribution of the probability of 
the A+ — > p + K~ + 7r + decay for the polarized A+ and the unpolarized proton are 
calculated in the effective quark model with chiral c7(3) x U(3) symmetry incorpo- 
rating Heavy Quark Effective theory (HQET) and Chiral perturbation theory at the 
quark level (CHPT)q. The application of the obtained result to the analysis of the 
polarization of the A+ produced in the processes of photo and hadroproduction is 
discussed. We draw the similarity between the measurements of the polarization of 
the A+ in the A+ —* p + K~ + ir + decay and the ^"-meson in the \x~ — > e~ + ^ e + ^ 
decay. 
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1 Introduction 

It is rather likely that in the reactions of photo and hadroproduction the charmed baryon 
A+ is produced polarized [1]. The restoration of the A+ polarization by means of the 
investigation of the decay products should clarify the mechanism of the charmed baryon 
production at high energies. 

The most favourable mode of the A+ decays to be detected experimentally is A+ — > 
p + K~ + 7r + . The experimental probability of this mode equals to -B(A+ — > pK~~-K + ) cxp = 
(5.0 ±1.3) % [2]. However, from the theoretical point of view this mode is the most difficult 
case due to the impossibility to factorize the baryonic and mesonic degrees of freedom for 
the computation of the matrix element of the transition [1]. 

The theoretical investigation of non-leptonic decays of charmed baryons without fac- 
torization of baryonic and mesonic degrees of freedom can be carried out in the effective 
quark model with chiral U(3) x U (3) symmetry incorporating Heavy Quark Effective The- 
ory (HQET) [3,4] and Chiral perturbation theory at the quark level (CHPT) g [5]. The 
effective quark model with chiral U(3) x U(3) symmetry is motivated by the low-energy 
effective QCD with a linearly rising interquark potential responsible for a quark confine- 
ment [6]. The application of this model to the description of the low-energy properties 
of charmed mesons: mass spectra [7], coupling constants [7-9], the form factors of the 
semileptonic decays [10] and the probabilities of the decays [7,11] gave the results agreeing 
good with experimental data. 

Recently [12] the effective quark model with chiral U(3) x U(3) symmetry has been 
extended by the inclusion of the low-lying baryon octet and decuplet coupled with the 
three-quark currents. Due to the dynamics of strong low-energy interactions caused by 
a linearly rising interquark potential there has been shown [6] that: (i) baryons are the 
three-quark states [13] and do not contain any bound diquark states, then (ii) the spinorial 
structure of the three-quark currents is defined as the products of the axial-vector diquark 
densities q C i(x)^qj(x) and a quark field qk(x) transforming under SU(3)f x SU(3) C group 
like (6/,3 c ) and (3y,3 c ) multiplets, respectively, where i,j and k are the colour indices 
running through % — 1, 2, 3 and q = u,d or s quark field. This agrees with the structure 
of the three-quark currents used for the investigation of the properties of baryons within 
QCD sum rules approach [14]. The fixed structure of the three-quark currents allows to 
describe all variety of low-energy interactions of baryon octet and decuplet in terms of 
the phenomenological coupling constant g-Q describing the coupling of the baryon octet 
and decuplet with the three-quark currents [12]: 

£intO*0 = ~^9b Bs(x) r]s(x) + g B B m (x) r] m (x) + h.c, (1.1) 

where B%(x) = (ifj p (x), . . .) and Biq(x) = (ipA++, ■ • •) are the fields of the baryon octet 
and decuplet, respectively, and r)s(x) = (—e ljk [u c i (x) / y^Uj(x)}'y l _ l 'y 5 dk(x) , . . .) and rjio(x) = 
(e tjk [u c i (x) / y' M Uj(x)]uk(x) , . . .) are the three-quark currents. The numerical value of pe, 
calculated in terms of the coupling constant svnn — 13.4 of the 7rNN interaction, has 
been found equal to #b = 1-34 x 10~ 4 MeV [12]. The coupling constants svna and <? 7 na 
of the 7rNA and 7NA interactions relative to the coupling constant svnn and the a^- 
term of the low-energy 7rN-scattering have been calculated in good agreement with the 
experimental data and other phenomenological approaches based on QCD [12]. 
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In this paper we apply the effective quark model with chiral U(3) x £7(3) symmetry 
to the calculation of the proton energy spectrum of the A+ — * p + K~ + 7T + decay of the 
polarized A+ and the analysis of the restoration of the A+ polarization in the dependence of 
the energies and momenta of the decay products. For the analysis of the A+ — > p+K~-\-ir + 
decay at the quark level we assume that the A+ is a three-quark state coupled with the 
three-quark current r] A +(x) = —e^ k [u c i(x)'y fl dj(x)}'j II '-f 5 Ck{x) defined as the product of the 
axial-vector light-diquark density u c i{x)^dj{x) and the c-quark field Cfc(x). 

The paper is organized as follows. In Sect. 2 we discuss the effective low-energy La- 
grangian describing non-leptonic decays of charmed hadrons and reduce the calculation 
of the amplitude of the A+ — > p + K~ + ir + decay to the calculation of the matrix element 
of the low-energy transition A+ — > p + K~ . In Sect. 3 we calculate the matrix element 
of the low-energy transition A+ — > p + K~ . In Sect. 4 we calculate the proton energy 
spectrum of the A+ — > p + K~ + n + decay for the polarized A+ and the unpolarized 
proton. In Sect. 5 we calculate the probability of the A+ — > p + K° decay relative to the 
probability of the A+ — > p + K~ + 7r + decay. The theoretical result on the ratio of the 
probabilities agrees well with the experimental data. In the Conclusion we discuss the 
obtained results. In the Appendix we calculate the momentum integrals describing the 
the matrix elements of the low-energy transitions A+ — > p + K~ and A+ — > p + K°. 

2 Effective Lagrangian for non— leptonic decays of charmed 
hadrons 

The effective low-energy Lagrangian responsible for non-leptonic decays of charmed 
hadrons reads [15] 



where Gf = 1-166 x 10~ 5 GeV -2 is the Fermi weak constant, V* s and V u d are the elements 
of the CKM-mixing matrix, Ci(fi) (i = 1,2) are the Wilson coefficients caused by the 
strong quark-gluon interactions at scales p > fi (short-distance contributions), where \i 
is a normalization scale. In the absence of quark-gluon interactions the coefficients C\{p) 
and C 2 (/i) do not depend on /i and amount to C\ — 1 and C 2 = 0. In (CHPT) g we 
should identify fi with the scale of spontaneous breaking of chiral symmetry (SB^S) A x = 
940 MeV [5-11], i.e. fj, = A x = 940 MeV. Therefore, below we would deal with Ci(A x ) 
and C 2 (A X ). The contribution of strong low-energy interactions at scales p < fi = A x 
(long-distance contributions) is described by (CHPT) g in terms of constituent quark loop 
diagrams, where the momenta of virtual quarks are restricted from above by the SB%S 
scale A x [5-12]. 

The structure of the term proportional to C 2 (A X ) can be reduced to the structure of 
the first one by means of the Fierz transformation [15]. The resultant coefficient of the first 
term would contain C 2 (A X ) in the form C 2 (A X )/A^, where = 3 is the number of quark 
colours. Thus, the effective low-energy Lagrangian responsible for the A+ — > p + K~ + 7r + 



£ e ff(x) 



G f 

V2 



V c * s [CM [s{x) y(l - 7 5 ) c(x)\ [u(x) - 7 5 ) d(x)} 

+CM [u(x) - 7 5 ) c(x)} [s(x) 7/i (l - 7 5 ) d(x)]}, (2.1) 
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decay can be taken in the form 

C cS (x) = V c * s V ud C l (A x ) [s(x)l,(l -l 5 )c(x)} i&ru(x)Y(l-l 5 )d(x)}, (2.2) 

where Ci(A x ) = Ci(A x ) + C2(A x )/3. The amplitude of the A+ — > p + K~ + n + decay is 
then defined 

where Ei (i = A^,p, K~, tt + ) are the energies of the A+, the proton and the K~ and 7r + 
mesons, respectively, and V is the normalization volume. 

In (CHPT) g at the tree-meson approximation we can factorize the pionic degrees of 
freedom and represent the amplitude Eq.(|2.3j) as follows 



2E A +V2E P V2E K -V2E V+ V 

= V* K^i(A x ) < p(q)K-(q_)\\s(0) 7 „(1 - 7 5 ) c(0)|A+(Q) > 

x <7r+(g + )|M(0)7 /i (l-7 5 )^(0)|0>, (2.4) 

where the matrix element < ir + (q + )\u(0) 7^(1 — 7 5 ) d(0)|0 > can be expressed in terms of 
the leptonic constant of the 7r + -meson F n = 92.4 MeV [5]: 



2E n+ V < n + (q + )\u(0) 7 "(1 - 7 5 ) d(0)|0 >=%V2F 7T <£. (2.5) 

Substituting Eq.( [2.5|) in Eq.( [2.4| ) we arrive at the amplitude of the A+ — > p + K~ + n + 

decay 

M(AUQ)^p(q)K-(q-)7r + (q + )) 



I2E A +V2E P V2E K -V 
= —%G F V* s V ud F w C 1 (A X ) ql < p(q)K~(q_)\s(0) 7m (1 - 7 5 ) c(0)|A+(Q) >, (2.6) 

where the matrix element < p(q)K~ (g_)|s(0) 7^(1 — 7 5 ) c(0)|A+(Q) > describes the tran- 
sition A+ — > p + K~ induced by the current s(0)j^(l — 7 5 )c(0) and defined by strong 
low-energy interactions. 

3 Strong low— energy transition A+ —> p + K~ 

By applying the reduction technique we bring up the matrix element of the strong low- 
energy transition A+ — > p + K~ to the form 



2E At V2E p V2E K -V2E w+ V < p(q)K'(q_)\s(0) lfl (l - 7 & ) c(0)|A+(Q) >= 

lim % [ d 4 Xl d A x 2 d A x 3 e l< i ■ x ^ e l< i- ' x * e~ l ® ' x * uJq, a' ) 

Q^Ml + ,q 2 ^Mlq 2 _^M 2 K J 
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\vf—-M p j{n 2 + Mj [ ) < 0|T(^(x 1 )^-(x 2 )[ S -(0) 7 ^(l -7 5 )c(0)]^ A +(x 3 ))|0 > 

(-il a -^*-M At y At (Q,*), (3.1) 

where ip p (xi), ifK-{ x 2) and Tp A +{x 3 ) are the operators of the proton, the i^~-meson and 
the A+ interpolating fields, u p (q, a' ) and u A +(Q, a) are the Dirac bispinors of the proton 
and the A+, respectively. 

Following [6-11] in order to describe the r.h.s. of Eq. ( |3.1p at the quark level we suggest 
to use the equations of motion 



V 7 "^ ~ Mp ) = vf 7/N ^' 

(n 2 + M 2 K )^ K -(x 2 ) = ^u(x 2 )z 7 5 s(x 2 ), 

^At(^)[-ir^-M At ^^fj At (x 3 ). (3.2) 

Here and gc are the phenomenological coupling constants of the proton and the 
A+ coupled with three-quark currents ^(^1) = — £ ufc [^ c j(^i)7 M ^j(^i)]7/x7 5 ^fc(2a) and 
= ^ fc Ci(x3)7 At 7 5 [J j (x3)7^^(x 3 )] [11], respectively: 

45? 0*0 = ^| $p(x) Vn(x) + ^= V A +(x) ip A +{x) + h.c. (3.3) 

Then, i,j and k are colour indices and u c (x) = u(x) T C and C = —C T = — = —C^ 1 is 
a matrix of a charge conjugation, T is a transposition. The three-quark current r] A +(x) = 
—e l ^ k [u c i{x)^ dj{x)} r ) ^ Ck{x) coupled with the A+ is constructed by analogy with rj^(x) = 
—e l ^ k [u c i{x)^ >J, Uj{xi)\^ ll ^ b dk{x) due to the similar spinorial properties of the A+ and the 
proton. Then, M p = 938 MeV and M A + = 2285 MeV are the masses of the proton and 
theA+. 

The interaction of the iT~-meson with quarks is described by the coupling constant 
gxgq = y/Zm/F^ 

£${x) = ^ s{x) i^uix) <p K - (x) + . . . + h.c, (3.4) 

where m = 330 MeV is the constituent quark mass calculated in the chiral limit [5] . 
Substituting Eq. (|3.2[) in Eq.( p.l|) we obtain 



2E At V2E p V2E K -V2E n+ V < p(q)K~(q_)\s(0) Tm (1 - 7 5 ) c(0)|A+(Q) >= 

= gB gc i ™ f d^x^xz e i( i ■ x ± e i( i- ■ x * e~ % Q ' x * 
2 F n J 

u p (q,a') < 0|T(r ?N (x 1 )[M(x 2 )^ 5 s(x 2 )][ S -(0) 7 ' 1 (l -7 5 )c(0)]r ?A +(x 3 ))|0 > u Ai (Q,a), (3.5) 

where the external particles are kept on-mass shell, i.e. Q 2 = M 2 +, q 2 = M 2 and 
q 2 - = M 2 K . 



5 



By applying the formulae of quark conversion [5] (Ivanov) we can determine the vac- 
uum expectation value in Eq. ( |3.3| ) in terms of the constituent quark diagrams. In the 
momentum representation we get [5-11]: 

2E A+ V2E p V2E K -V2E 7r+ V < p(q)K-(qJ)\s(0) 7 „(1 - 7 5 ) c(0)|A+(Q) >= 



-Wb9c 



3m 


1 


2 


fd 4 h r 




IQtt 2 


J 


n 2 i J 



n 2 i 



7 



-7m(1-7 5 )" 



u p (q, a' ) 7a7 5 
1 



-r 



-r 



m — h\ m + k 2 m — q + k\ + hi 
7Z37 5 ' %+(<?, (3.6) 



m — q — q- + &i + k 2 M c — Q + ki+k 2 

where M c = 1860 MeV is the mass of the constituent c-quark [7-1 1] . In the HQET the 
r.h.s. of Eq.(|J) is given by [3,4,7-11]: 



2E A+ V2E p V2E K -V2E w+ V < p(q)K-(q_)\s(0) 7 „(1 - 7 5 ) c(0)|A+(Q) >= 



~Wb9c 



3m 


1 


2 


fd 4 h r 


~K 


167T 2 


J 


7C 2 i J 



TT 2 i 



"7 



P. 



-r 



i 



7 



m — q — g_ + k\ + k 2 



7m(1-7 5 ) 



1+0 



m — k\ m + k 2 m — q + k\ + k 2 
1 



f7/3 7 u A +(Q,a), 



(3.7) 



2 / [(fei + k 2 ) ■ v + z'0] 

where is 4-velocity of the A+ [3,4,7-11] normalized by v^v^ = 1. 

For the computation of the momentum integral we assume [12] that the proton is a 
very heavy state and its 4-momentum is much larger than other momenta in the integrand 
of Eq. Q3.7p . Keeping the leading terms in the large M p expansion we reduce the r.h.s. of 
Eq.(|3T7D to the form 



2E At V2E p V2E K -V2E n+ V < p(q)K- (q_)\s(Q) lp (l - 7 5 ) c(0)|A+(Q) >= 



■ 9b9c 3m 
' l M 2 F w 

l+v\ 



1 


2 


rd 4 h 


r d 4 k 2 


167T 2 


J 


TT 2 i J 


n 2 i 



l 



l 



-t r 

m — k\ m + k 2 



7 a 7 ,(l- 7 5 ) 



f7/3 7 5 u A +(Q,a) 



\ 2 J [(fci + k 2 ) ■ v + iO] 
The replacement of the constituent quark Green function 

1 Pi 

Ml 



(3.8) 



(3.9) 



m - q + h + k 2 

agrees with the heavy baryon [16,17] and HQET [3,4] approaches. Indeed, in accordance 
with [16,17] and HQET [3,4] we obtain 

1 m + q — k\ — k 2 



m — q + k\ + k 2 



(q-h- k 2 y 



iO 



m + q — k\ — k 2 



M 2 - 2(fci + fc 2 ) • g - m 2 + (A* + k 2 ) 2 + zO 



1 



m — k\ — k 2 



+ v 



M p 2(fci + fc 2 ) ■ v _rn^, (h + k 2 f ' 
M p M 2 + M 2 + ' U 



(3.10) 
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where we have set = M p [3,4,16,17]. In the case of m = M p and in the limit M p oo 
we arrive at the well-known expression for the Green function of a heavy baryon (or a 
heavy quark in HQET [3,4]) [16,17]: 



1 1 



m — q + k\ + &2 M p + k\ + k 2 — M p v 

fl + v\ 1 



(fci + k 2 ) -v + iO' 



(3.11) 



In our case m <C M p , therefore, in the limit M p — > oo [16,17] we arrive at the expression 
Eq.(H). 

The computation of the momentum integrals in Eq.( |3.8| ) is carried out in the Appendix. 
Thus, the matrix element < p(q)K~ (g_)|s(0) 7 M (1 — 7 5 ) c(0)|A+(Q) > is defined 



2£ A +V2£ p V2 J EW2£ 7r +V <p(g)A-(g_)|s(0) 7M (l-7 5 )c(0)|A+(Q) >= 

= ig^ss ~ ^ fip(g, ^ )7«7 5 ^7 a 7,(l - 7 5 ) (^p) 7/.7 5 % C+ (Q, a), (3.12) 

where ^nn is the coupling constant of the 7rNN interaction expressed in terms of the 
parameters of the model through the relation svnn — 9B(2iTi/3F n )(< qq > 2 /M%) [12]. 
After some algebra the matrix element Eq. ( |3.12| ) amounts to 



/2E At V2E p V2E K -V2E n+ V < p(g)tf-( g _)|s(0) Ta4 (1 - 7 5 ) c(0)|A+(Q) >= 

= iftrNN t — ^' ) t 2 M 1 ~ 7 5 ) + 7m(1 + 7 5 )] u A + {Q, a) = 

= ig^nt — ^u p (q,a') (1 - 7 5 ) (2 + 7^) u A +{Q, a), (3.13) 
5 g-Q mr c 

where we used the Dirac equation of motion vu A +(Q, a) = u A +(Q, a) for the free A+. 



4 Proton energy spectrum of the A+ — ► p + K + 7r + 
decay 

The amplitude of the A+ — > p + A' - + 7r + decay is given by 

Al(A+(g) ^p(g)A-(g_)7r + (g + )) = G^Kd^AJ 

,/( ' ^ * 1 u p (q, a') (I- 7 5 ) (2Q ■ q + + M A +q + ) u A +(Q, a). (4.1) 



4 ^ttNN 

X 5% 



<7b "i 2 



The partial width of the A+ —>■ p + K + 7r + decay determined in the rest frame of the 
At reads 



dT(At ^pK-7r + ) = ^ f -\M(At(Q)^p(q)K-(q^(q + W 

(2 ^^. w _ 9 _ 9 __ 9+) ^_||__gt_, (4 ,) 
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We define the quantity \A4(A+(Q) — > p(q)K^(q_)ii + (q + ))\ 2 for the polarized A+ and the 
unpolarized proton: 

"1 2 

4 ^ttNN gC 



\M(A+(Q)^p(q)K-(q„)7r+(q + ))\* = V^A, 



5 M A + g B m 2 



x - tr{(Q + M A +)(1 + 7 5 ^a+)(2Q ■ g+ + M A +g + )(l + 7 5 )(g + M p )(l - T 
x (2Q-g + + M A +g+)}, 



(4.3) 



where oo^ + is the space-like unit 4-vector (u;1 



-1) orthogonal to the 4-momentum of 



the A+ iyj A + ■ Q = 0) and related to the direction of its spin. In the rest frame of the A+ 
we have: u>^ + = (0, oj a +) such as lj 2 + = 1. 

Neglecting the terms proportional to M 2 and M\ the result of the calculation of the 
trace reads 

itr{...} = 16(g-g + ) 2 (Q-g) + 24M^(Q.g + )(g. g+ )-32M A +(g-g + ) 2 K r g) 

+ 16M A +(Q ■ q + )(Q ■ q)(u At ■ q+) + 8M|+(g • q+)(u At ■ g + ). (4.4) 
For the derivation of the proton energy spectrum it is convenient to use the formula 

d 3 g_ d 3 q + 



J q a + qi5 {4 \Q-q-q--qA 



2Ej(- 2E n -t 



^ (~(Q ~ q) 2 if* + 4 (Q - q) a (Q - q?) J -q-q- 



d 3 q_ d?q + 
2Ex- 2E n + 



71 

24 



{Q-q) 2 9 aP + ±{Q-q) a {Q-qy 



(4.5) 



which is valid when the contributions proportional to M 2 and M\ are neglected. Using 
Eq.flPQ we get 

d 3 q_ d 3 q + 



J ± t r{...}5^(Q-q-q_-q, 



TV 

3 



2Ej{- 2E W + 

(15M A +(Q ■ q) - 18M 2 +(Q ■ q) 2 + 8(Q ■ qf + 7M 2 +M 2 (Q • q) - 12M A +M 2 ) 
A + • g)M A +(13M Ac+ - UM 2 At (Q • g) + 8(Q • g) 3 - 7.1/;. M 2 p ] 



—Ml + x 
2 A - 



' 3 

1 x 

5 



15 



x 2 + -£ 2 - 
60 s 



2^ 

5 x 



+ ^ 



13 



n„) — \1 r 1 xH x £ 



X 



13 



13 



52 



(4.6) 



where the final expression is taken in the rest frame of the A+, x = 2E p /M A + is the scaled 
proton energy, £ = 2M p /M A + and n p = q/\q\. The scaled proton energy x ranges the 
region £ < x < 1 + £ 2 /4. 

The proton energy spectrum of the A+ — > p + If - + 7r + decay in the rest frame of the 
A+ is determined: 



|g f \/;k^i(a x )P 



A'ttNN 



4gcKA> 
5 5> B ra 2 



5M I+ 

5127T 3 



X 1 



■ X + 



15 



60 S 



2_e 

5 x 



1 + a(x) (u A + ■ n p ) 



x \ x 



47T ' 



(4.7) 



where dD.fi is the solid angle of the unit vector n p = q/\q\ and a(x), the parameter of 
the asymmetry related to the polarization of the A+, is given by 

, 7 2 2 7 2 

n / c2 1 x H x 1 



3 2 2 7 ~ 2£ 

-H x 2 H £ 

5 15 60 s 5 x 



(4i 



1 



In order to apply Eq. (|4.7| ) to the analysis of the polarization of the A+ in the processes 
of photo and hadroproduction we have to remove an uncertainty related to the arbitrary 
coupling constant g 2 ^. For this aim we suggest to normalize the proton energy spectrum 
to the partial width of the mode A+ — > p + K~ + n + and replace the coupling constant 
by the experimental value of the probability. Integrating Eq. ( |4.7|) over all variables we 
obtain the partial width of the mode: 



4 gc F W A, 



5 g B m 2 



x 



5M|+ 

5127T 3 



x /(0- (4-9) 



The function /(£) is determined by the integral 

!+« 2 /4 , 



/(0 



7 



2f 



1 - - x + — x^ + — C ~ 7 — \x Jx 2 - £ 2 dx = 0.065 



15 



60 



5 x 



(4.10) 



The numerical value has been obtained at M A + = 2285 MeV and M p = 938 MeV. 
The proton energy spectrum of the A+ — > p + X~ + 7r + decay is then given by 



dB(A+ ^pK-n- 
B(A+ ^ pK-n+ 



7_ 
60 



15.40 x 1 



1 + a(x) (tu A + ■ n p ) x J x 2 — S, 2 dx 



2f 
5 x / 



47T 



By using the experimental value of the probability £>(A+ — > pK 7r + ) cxp 



(4.11) 
(0.050±0.013) 



[2] we obtain 

dB(A+ ^pR-^ 



(0.77 ±0.20) x (I -^x + ^x 2 + ^i 2 



2 

15 



7 

60 



5 x 



1 ± a(x) (cu A + ■ n t 



X \ X 



£ 2 dx ^. 

47T 



(4.12) 



Integrating over the energy of the proton we derive the angular distribution of the prob- 
ability of the A+ — > p + + 7r + decay: 



ciB(A+ -> p R-7V+) 0.050 ± 0.013 



47T 



1 + 0.77 (cu, 



(4.13) 



Thus, the proton energy spectrum Eq. (|4.12|) and the angular distribution of the probabil- 
ity Eq. (f4.13 ) of the A+ — ► p + K~ + n + decay do not contain arbitrary parameters and, 
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therefore, can be applied to the analysis of the polarization of the A+ in the processes of 
photo and hadroproduction. 

The formulae Eq. (|4.12|) and Eq. (|4.13| ) define the polarization of the A+ relative to the 
momentum of the proton. If the spin of the A+ is parallel to the momentum of the proton, 
the right-handed (R) polarization, the scalar product uj a + ■ n p amounts to uj a + ■ n p = cos 
The angular distribution of the probability reads 

dB(A+ ^pK-7T + ) {R) _ 0.050 ±0.013 



dVtf 



An 



x (1 + 0.77 cos i 



(4.14) 



In turn, for the left-handed (L) polarization of the A+, the spin of the A+ is anti-parallel 
to the momentum of the proton, the scalar product reads (uj a + ■ n p ) = — cost? and the 
angular distribution becomes equal 



dB(A+ ^ P K-n+) {L) _ 0.050 ±0.013 



dVtf 



An 



0.77 cos( 



(4.15) 



For the right- and left-handed polarizations of the A+ the proton energy spectrum is 
given by 



dB(A.+ ^pK~n + ) 



(R,L) 



(0.77±0.20)x i_ ,+ ,»+ {.- L 



X 



1 ± a(x) cosfi 



x \l x 2 — £ 2 dx — — . 
v An 



(4.16) 



The formulae Eq.( 4.12j ) and Eq.( 4.13| ) resemble the electron energy spectrum and the 
angular distribution of the probability of the /3-decay of the /i~-meson, i.e. \i~ — > e~ + 
v e + Vfi- Therefore, the procedure of the restoration of the polarization of the A+ in the 
A+ — > p + K~n + decay is in complete analogy to the procedure of the measurement of 
the polarization of the /!~-meson in the fi~ — > e~ ± v e + decay [2] . 



5 Probability of the A+ — > p + K decay 

Most modes of the A+ decays are measured relative to the mode A+ — >• p + K~ + n + [2]. 
For the theoretical description of the A+ — > p + K~ + n + decay at the quark level we have 
introduced the phenomenological low-energy interaction of the A+ with the three-quark 
current f] A +(x) containing an arbitrary phenomenological coupling constant gc Eq. ( |3.3|) . 
The spinorial structure of the three-quark current r) A +(x) = —e^ k [u c i(x)-y fl dj(x)]-y^'y 5 Ck(x) 
defined as the product of the axial-vector light diquark density [u c i(x)'j fl dj(x)] transform- 
ing under SU(3)f x SU(3) C group like (6^,3^,) and the c-quark field Cfc(x) is caused by 
the dynamics of the quark confinement given by a linearly rising interquark potential 
[6,11]. In order to verify the validity of the approach applied to the computation of 
the proton energy spectrum and the angular distribution of the probability of the mode 
A+ — > p + K~ + n + it is not sufficient to be restricted by the consideration only this mode. 
For the confirmation of the result obtained for the mode A+ — > p + K~ + n + one needs the 
computation of the probabilities of other modes relative to the probability of the main 
mode A+ — > p + K~ + n + . In the ratio the coupling constant gc cancels itself and the 
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theoretical result turns out to be dependent on the Wilson coefficients, determined by 
the short-distance quark-gluon interactions, and the long-distance dynamics, described 
by the effective quark model with chiral U(3) x U(3) symmetry motivated by QCD with 
a linearly rising confinement potential. The agreement between the experimental data 
and the theoretical predictions for the ratios should testify both the self-consistency of 
the approach and the consistency of it with a short-distance QCD. Below we obtain the 
evidence of the self-consistency and the consistency of the approach by example of the 
calculation of the decay mode A+ — > p + K°. We have chosen this mode due to the 
following reasons. First, it contains the proton in the final state, and, second, the com- 
putation of the matrix element of this mode is completely different to the computation 
of the matrix element of the mode A+ —>■ p + K~ + 7r + . Indeed, the mode A+ —>■ p + K° 
unlike the mode A+ — > p + K~ +7r + admits the factorization of the baryonic and mesonic 
degrees of freedom. 

The effective low-energy Lagrangian responsible for the decay A+ — > p + K° can be 
obtained from the effective Lagrangian Eq.( [2.1| ) at fi = A x : 



C cS (x) = V c * s V ud C 2 (A X ) [u(x) 7,(1 - 7 5 ) <x)} [S{x)-f(l ~ 7 5 ) d(x)), (5.1) 
where C 2 (A X ) = C 2 (A X ) + C^/Z. 



The amplitude of the A+ —>■ p + K° decay can be defined by analogy with Eq.( p.6|) : 



where Fjc = 113MeV [2] is the leptonic constant of the ii'-mesons. 

To the computation of the matrix element < p(q)\u(0) 7,(1 — 7 5 ) c(0)|A+(Q) > we 
apply the reduction technique. By using the equations of motion Eq. (|2.2| ) we arrive at 
the expression 



x u p {q,a') < 0|T(77n(xi) [m(0) 7 m (1 - 7 5 )c(0)]r /Ac+ (x 2 ))|0 > u A +{Q,a), (5.3) 



M(Aj(Q)^p(q)K (q )) _ 
y /2E A +V2E p V 

= -%G F V c * s V ud C 2 (A x ) F K q» < p(q)\u(0) 7 ,(1 - 7 5 ) c(0)|A+(Q) >, 



(5.2) 



pE At V2E p V < p(q)\u(0) 7 „(1 - 7 5 ) c(0)\K(Q) > 




where the A+ and the proton are kept on-mass shell, i.e. Q 2 = M\ + and q 2 = M 2 In 
terms of the constituent quark diagrams represented by the momentum integrals the r.h.s. 
of Eq. ( p.3[ ) defined in HQET reads 



y/2E A +V2E p V <p(q)\u(0)l^-l 5 )c(0)\At(Q) >= 




(5.4) 
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Keeping the leading terms in the large M p expansion we reduce the r.h.s. of Eq.( |5.4|) to 
the form 



2E A +V 2E P V <p(q)\u(0)l^-r)c(0)\K(Q) >= 

1 2 -'41. „ M I 



9b9c 

M 2 
p 

'1 + v 



16vr s 



<rk x r drk 2 _ , , . e 1 
)7a7 



7Tn J 7T"? 



-y 



■P. 



m — k\ m + k 2 



7 a g 7 M(l-7 5 ) 



(5.5) 



The integrals over k\ and k 2 have been calculated in the Appendix. Using Eq.(A.7) and 
making some algebraic transformations with the Dirac matrices we get 



l2E Kt V2E p V < p(q)\u(0) 7/ ,(l - 7 5 ) c(0)|A+(Q) >= - ^ 
x u p (q, a')(2v-f fl q + vgj^l - 7 5 ) %+(<?, 
The amplitude of the — > p + decay is defined 



NN 



4 £c FA 



5 (7b 



(5.6) 



A4(A+(Q) ^p(q)K°(q )) = iG F V^V mi C 2 {A x )^- 

m 

xu p (q,a')(A + B^ 5 )u A +(Q,a), 
where the constants A and B are given by 

A = 1 
B = 1 



4 g c F T A, 
SVnn 7 j- 

5 pb m 



(5.7) 



M p 


- 2 


M 2 
p 


M| 


M At 




M h 


Kt 


M p 


- 2 


M 2 

p 


M| 


M At 




Kt 


Kt 



1.03, 
0.21. 



(5i 



The numerical values are obtained for M A + = 2285 MeV, M p = 938 MeV and M R0 
498 MeV [2]. The partial width of the A+ 4 p + K° decay reads 



r(A+ - pK°) = \Gp V* V ud C 2 (A x )\ 2 Lnn \ — ^y 

5 (7b 

x tr{(Q + M A+ )(A - B^)(q + M P )(A + 5 7 5 )} 



2 F% 



m? 327T 



\ 



1 - 



'M p - M R0 \ 



1 



' M P + M R0 

m7+ 



(5.9) 



When computing the trace over Dirac matrices we obtain the partial width in the form 

4 g c F W A X 



T(At ^ P K°) = \G F V; s V ud C 2 (A y 



9nNN 



32vr 2 F, 



x 



2 t?2 
K 



5 m 2 



-4 < A 1 



M\ 



K 



5 g B rn 2 
+ B 2 



' ^Kt 



5127T 3 
\ 2 
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X 



U 



1 - 



M p — M 



2i 



K° 



M 



1 - 



M v + M, 



K° 



11.72 x \G F V; s V ud C 2 (k, 



9nNN 



4 gc F v A y 



5 gs rn? 



2 r 



5M 5 A+ 

Il c 

5127T 3 



(5.10) 



Now we can define the partial width of the A+ — > p + K° decay with respect to the partial 



width of the A+ — > p + K + 7r + decay. By using Eqs. ( |4.9| ) and ( [4.1 0| ) we get 



R 



th 



r(A c + -»pjf°) 
T(A+ -> pK n 



180.31 x 



gf(A> 



(5.11) 



The numerical factor in front of the ratio of the squared Wilson coefficients is completely 
due the low — energy dynamics of our model induced by a linearly rising confinement 
potential. In order to verify the consistency of this dynamics with a short-distance QCD 
we should substitute in Eq. ( |5.11|) the numerical values of the Wilson coefficients. 



Following Buras et al. [15] we obtain Ci(A x ) = 1.24 and C 2 (A y 



-0.47. These 



numerical values agree well with the numerical values of the Wilson coefficients calculated 
at the normalization scale p ~ 1.5 GeV [15]: C\(p) ~ 1.21 and C 2 (p) — —0.42. 



The ratio R t h calculated at Ci(A x ) = 1.24 and C 2 (A y 



-0.47 amounts to 



R 



th 



r(A c + -> P jr°) 

T(A+ -> pK-n 



0.50. 



(5.12) 



The theoretical result agrees well with the experimental value averaged over all exper- 
imental data [2]: -R ex p = (0.49 ± 0.07). This agreement is non-trivial and testifies not 
only the consistency of the model with short-distance QCD but the self-consistency of 
the effective quark model with chiral U(3) x U(3) symmetry incorporating HQET and 
(CHPT)^. The former is due to the distinction between the computations of the matrix 
elements of the modes A+ —>■ p + K° and A^ —>■ p + K~ + 7r + . Indeed, the computation 
of the matrix element of the mode A+ — > p + K° admits the factorization of the baryonic 
and mesonic degrees of freedom, whereas for the computation of the matrix element of 
the mode A+ — > p + K~ + 7r + such a factorization is not feasible. 

Our theoretical result for ratio Eq. (|5.12|) confirms too the validity of our prediction 
for the proton energy spectrum and the angular distribution of the probability of the 
A+ — > p + K~ + 7r + decay given by Eq. (|4.12|) and Eq. (|4.13|) , respectively. 



6 Conclusion 

The main result of the paper is in the prediction of the proton energy spectrum and the 
angular distribution of the probability of the mode A+ — >• p + K~ + n + of the A+ decays. 
This mode is the most favourable for the measurement as it contains the proton and 
the charged mesons. However, from the theoretical point of view this mode is the most 
difficult for the computation due to the impossibility to factorize baryonic and mesonic 
degrees of freedom. 

To the computation of the matrix element of the A+ — > p + K~ + 7r + decay we have 
applied the effective quark model with chiral U(3) x £7(3) symmetry incorporating HQET 
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and (CHPT)g. This model is motivated by the effective low-energy QCD with a linearly 
rising confinement potential. Due to the dynamics of strong low-energy interactions 
caused by a linearly rising confinement potential the spinorial structure of the three- 
quark currents coupled to the baryons is fixed unambiguously. The effective low-energy 
interactions of the low-lying baryon octet and charmed baryons coupled to the three- 
quark currents can be described in terms of two phenomenological coupling constants g-Q 
and gc, respectively. These constants enter multiplicatively to the matrix elements of the 
strong low-energy transitions of baryons. 

In the case of the non-leptonic decays of the A+ the multiplicative character of the 
constants g-Q and gc allows to replace the product of these constants by the experimental 
value of the probability of the mode A+ — > p + K~ + 7r + . This defines any mode of the A+ 
decays relative to the mode A+ — > p+K~+n + . As regards the proton energy spectrum and 
the angular distribution, the resultant expressions do not contain any arbitrary parameters 
and can be applied to the analysis of the polarization of the A+ in the processes of photo 
and hadroproduction. We have considered the simplest case, maybe most favourable 
from the experimental point of view, when the A+ is polarized while the proton of the 
decay is unpolarized. This means that for the restoration of the polarization the A+ one 
should follow only the geometry of the momenta of the protons of the decay but not their 
polarizations. In this case there is an obvious similarity between measurements of the 
polarization of the A+ in the A+ — > p + K~ + ir + decay and the /z~-meson in the /3-decay 
(J>~ -> e~ + v c + Vp. 

For the confirmation of the validity of our prediction for the proton energy spectrum 
and the angular distribution of the probability of the mode A+ — > p + K~ + n + , we have 
computed the probability of the mode A+ — > p + K° relative to the probability of the 
mode A+ — * p + K~ + 7r + . Our prediction for the ratio of the probabilities R t h = B(Aq — > 
pK )/ B(A+ — > pK~ir + ) = 0.50 agrees well with the experimental value averaged over all 
experimental data -R ex p — (0.49 ± 0.07). This agreement is not trivial and confirms not 
only the self-consistency of our approach but the consistency of it with a short-distance 
QCD, since the computation of the matrix element of the mode A+ — > p + K° differs from 
the computation of the matrix element of the mode A^ — > p + K~ + 7r + . Indeed, if for 
the computation of the matrix element of the mode A+ — > p + K° one can factorize the 
baryonic and mesonic degrees of freedom, whereas in the case of the computation of the 
matrix element of the mode A+ — > p + K~ + 7r + such a factorization is not feasible. 
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Appendix. Computation of the momentum integrals 

We perform the integration over fa and k 2 of the momentum integral of Eq. ( ^.6[ ). For this 
aim we consider the integral 



J(v) 



d 4 h fd 4 k 2 1 1 



1 



ir 2 i J ttH m - fa m + k 2 [(fa + k 2 ) ■ v + iO] ' 



(a.i; 



where T is a Dirac matrix. The nontrivial contribution comes from the components of 
the k 1 ^ and k 2 4-vectors parallel to 4-vector v M . This gives 



J(v) 



d A ki r d A k 2 m + vfa ■ v m — vk 2 ■ v 1 

~nH J ~nH [m 2 - k\ - iO] [m 2 - k\ - io] [(fa + k 2 ) ■ v + iO] 



(A.2) 



Now it is convenient to make the Wick rotation and pass to Euclidean momentum space 
[9]: 



rrl \ A ■ f dk E lk 3 E1 f dfli /■ 

J(v) =A l j J^ 2 -(m + l v k E1 cos X i)r y 



dk E2 k% 2 
m 2 + fcl 



£2 



<if2 2 m — ivk E2 cos X2 
2ir 2 k El cosxi + k E2 cos\2 



(A.3) 



where dQj = 47r sin 2 Xi^Xi (i = 1, 2) are the solid angles in Euclidean spaces of the mo- 



menta k El and k E2 , respectively, k Ei = ykli + k?(i = 1,2). Then we have used the 

relation y^v% = —i [8]. 

In order to disconnect integrations over k E \ and /c£2 we suggest to use the following 
integral representation 



oo oo 







- (m - iv k E2 cos X 2 ) e l t km COS X2 ■ 



dk E2 k E2 f dVL 2 



m 2 + k 2 E2 J 2tt 2 



(A.4) 



Integrating out Xi an d X2 we get 

dk E \k E1 



00 00 



J(v) = 16 J dt J 







dk E2 k 



m 2 + k\ 



m 



El 



Ji(k m t) 



k E1 t 



~ lV dt 



Ji(k E1 t) 



k E1 t 



E2 



m 2 + k\ 



m 



E2 



Ji(k E2 t) 
k E 2t 



d 



Ji(k E2 t) 
k E 2t 



(A.5) 



where J\(k E {t) (i = 1, 2) is the Bessel function. Now we can perform the integration over 
kEi (i = 1,2): 



00 , 



K x (mt) 



, d 



K^mt) 



>r< m 



Ki(mt) 



, d 



Kx(mt) 



(A.6) 
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where Ki(mt) is the McDonald function. The integrals over t are divergent and should 
be regularized. We suggest to use the cut-off regularization restricting the region of the 
integration over t from below as t > 1/A X . Keeping leading divergent contributions [5-12] 
it is convenient to represent the r.h.s. of Eq.(A.6) in the following form: 



J(v) = 5 



16tt 2 



— | < qq> 2 vTv, 



(A.7) 



where < qq > is the quark condensate defined in terms of the SB^S scale and the con- 
stituent quark mass as follows [5-12] 



< qq > = 



N 



fd 4 k 



16tt 2 



m 



Nm 
"4^ 



A 2 

x 



The numerical value is calculated at N 



m 



-mHn[\ + -^\ =-(253MeV) 3 . 
\ m 2 J 

(A.8) 

330 MeV and A v = 940 MeV. As has 



been shown in Ref. [5] the quark condensate value < qq >= —(253 MeV) 3 describes with 
an accuracy better than 5% the mass spectrum of low-lying pseudoscalar mesons for the 
current quark masses m 0u = 4 MeV, m od = 7 MeV and m 0s = 135 MeV quoted by QCD 
[18]. 
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